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Abstract. Molecular clouds act as primary targets for cosmic-ray interactions and are expected 
to shine in gamma-rays as a by-product of these interactions. Indeed several detected gamma- 
ray sources both in HE and VHE gamma-rays (HE: 100 MeV < E < 100 GeV; VHE: E > 100 
GeV) have been directly or indirectly associated with molecular clouds. Information on the local 
diffusion coefficient and the local cosmic-ray population can be deduced from the observed gamma- 
ray signals. In this work we concentrate on the capability of the forthcoming Cherenkov Telescope 
Array Observatory (CTA) to provide such measurements. We investigate the expected emission from 
clouds hosting an accelerator, exploring the parameter space for different modes of acceleration, age 
of the source, cloud density profile, and cosmic ray diffusion coefficient. We present some of the 
most interesting cases for CTA regarding this science topic. The simulated gamma-ray fluxes depend 
strongly on the input parameters. In some cases, from CTA data it will be possible to constrain both 
the properties of the accelerator and the propagation mode of cosmic rays in the cloud. 
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INTRODUCTION 

The physical scenario Emission in HE-VHE y-rays is expected in spatial coincidence 
with molecular clouds, resulting from the hadronic interaction between cosmic-ray (CR) 
particles and the dense material in the cloud acting as a target. Indeed, some MCs have 
been detected in y-rays in both the GeV and TeV domain [see, e.g., 1, 2, 3]. Moreover, 
it has been suggested that some of the still unidentified y-ray sources might also be 
MCs illuminated by CRs that escaped from an accelerator located inside the cloud 
or in its proximity [4, 5, 6, 7]. In such cases the modeling of the emission involves 
the parametrization of the diffusion of charged particles. The study of this emission 
is extremely useful in unveiling the physics of CR sources. A more detailed study 
following the lines of the present contribution can be found in [8]. 

If a power-law energy spectrum (J p (E p ) = KE p y ) is assumed for the intensity of pri- 
mary CRs, the resulting y-ray spectrum due to hadronic interactions would also follow a 
power-law spectrum (F(E) E~ r ). The y-ray spectrum reproduces the spectrum of the 
parent particles, at high energies. Although, if we consider energy-dependent diffusion 
coefficient, the CR spectrum may differ from a simple power-law near the acceleration 



site. Here, we assumed that the source is point-like and located at the origin of the coor- 
dinate system. The acceleration and diffusion processes are computed following the ap- 
proach of [5], with the appropriate solution of the diffusion-loss equation. The diffusion 
coefficient is assumed to depend on the CR energy only, as: D(E) = D\q(E / lOGeV) 5 
cm 2 s _1 . We investigate the case of an accelerator at the center of a molecular cloud. 

The resulting flux in y-rays is mainly dependent on: diffusion coefficient and its 
energy dependence (Dio, <5); age of the accelerator; type of accelerator (impulsive/ 
continuous); spectrum of injection (y); fraction of energy in input (total energy in form of 
cosmic rays W p = r\ 10 50 erg). An impulsive source of particles corresponds to the case 
when the bulk of relativistic cosmic-rays are released during times much smaller than the 
age of the accelerator itself. When the timescales are comparable, the source is referred 
to as a continuous injector. The total energy input is taken as W p in the impulsive case, 
while the energy injection rate in the continuous case is of L p = 10 37 erg s _1 , resulting 
in the same total input for accelerators with r\ = 1 and with an age of a few hundreds 
of thousand years. In this scenario, a great variety of y-ray spectra is expected. This can 
produce a variety of different GeV-TeV connections, some of which could explain the 
observed phenomenology [see, e.g., 11]. 

CTA The Cherenkov Telescope Array (CTA) is an international project for the devel- 
opment of the next generation ground-based very high energy (VHE) gamma-ray instru- 
ment. CTA will significantly advance with respect to the present generation IACTs: it 
will feature an order of magnitude improvement in sensitivity at the core energy range of 
1 TeV, improve in its angular and energy resolution, and provide wider energy coverage, 
see [12]. Indeed, the array is expected to have an unprecedented sensitivity down to ~50 
GeV and above ~50 TeV, establishing a strong link with the satellite-based operations 
at low energies, namely the Large Area Telescope on board the Fermi satellite, see [13] 
and water Cherenkov experiments at the highest energies [e.g, HAWC, see 14]. Both a 
southern and a northern hemisphere observatory are foreseen. 

Detailed response functions for a proposed CTA array can be found in [15] . We use the 
procedure detailed there to calculate the spectral points and profiles from CTA simulated 
observations. 



Spectral features A constraint on the diffusion coefficient can come from the identifica- 
tion of a break in the y-ray spectrum integrated from the entire cloud region. The break 
can be related to the minimum energy that can diffuse in the entire cloud over a timescale 
comparable to the age of the accelerator. Equating the diffusion timescale, with the age 
of the accelerator, one obtains the constraints shown in Fig. 1, left, corresponding to: 



For scenarios with galactic average diffusion (Dio = 10 cm s ) and energy depen- 
dence parameter in the range 8 = [0.3. .0.6], the corresponding break in y-ray emission 
will always be at energies below the CTA energy acceptance. This is accurate for im- 
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pulsive accelerators. In the continuous acceleration case, new injections of high energy 
particles will smooth the effect of a break. At energies higher than the break given by 
Eq. 1, the particle spectrum will follow a power-law form composed by the slope of the 
injection spectrum and the energy dependence of the diffusion coefficient (i.e y + <5). 
Therefore the y-ray emission will also show a power-law behavior, reducing the capa- 
bility of constraining the parameter space from y-ray data. 

Let us assume that a molecular cloud with measured mass and distance is detected 
at TeV energies. Let us further assume that from multiwavelength observations we 
identified a possible accelerator of CRs responsible for the gamma ray emission and 
that an estimate of the age of that accelerator is known. In Fig. 1 we show the simulated 
spectra for the gamma ray emission for such an accelerator, which is assumed to be 
impulsive and with an age of 10 4 years. The molecular cloud is assumed to have a mass 
of M5 = 1O 5 M and a radius of 20 pc (hence with an average density of «h = 130 cm~ 3 ) 
located at a distance of d =1 kpc. The other parameters are varied on a discretized grid 
and the corresponding observed spectrum is simulated from CTA responses. The best 
fitting model belonging to the parameter space grid is found along with the models with 
a X 2 /dof not more different than unity from the best fit model. The right panel of Fig. 
1 shows the case of D m = 10 26 cm 2 s _1 , 8 = 0.4, y = 2.3, 77 = 1/3. Thanks to the high 
flux reached in this case and the presence of a break, the spectrum is reconstructed easily 
to the intrinsic parameters, with a break at E r ~ 2 TeV and slopes T\ = 2.3 and T2 = 2.7 
below and above the break, respectively. 
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FIGURE 1. Left Boundaries for £break > 70 GeV in y-ray (left from the boundary) for different ages of 
the accelerator (10 3 ' 4 ' 5 years are represented by the tick lines, solid, dashed, dotted curves, respectively). 
The thin solid line represents the boundary for an age of 10 4 years and £break > 1 TeV. The point 
represents the parameters of the spectrum given in the right panel. Right CTA expected performances 
on the reconstruction of the intrinsic model for D\q = 10 26 cm 2 s _1 , 8 = 0.4, y = 2.3, 7] = 1/3. The 
points are the expected spectral points from 50 hours of CTA observation time. The intrinsic spectrum 
is the only accepted model (see text). The expected 5 year point source sensitivity for Fermi/LAT is also 
given. This is calculated from the 1 year sensitivity in [13], linearly scaled with time at high energies 
(> 10 GeV). 

Morphology The different scenarios can also be disentangled by investigating the 
morphology and extension of the emission region. The shape of those profiles depends 
on the parameters of the simulated scenario, an example of which is shown in Fig. 2. 
These profiles are easily distinguishable from each other. Extensions of the y-ray emis- 
sion depend on the parameters studied here, with some general trends. Older sources 
are always more extended than younger sources, as the lower energy particles will have 
diffused further from the center of the cloud and thus from the accelerator. Emission due 
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FIGURE 2. Profiles of the photon count, normalized to the respective maximum. Here we show the 
example of impulsive (left) and continuous (right) acceleration, Diq = 10 26 cm 2 s _1 , 8 = 0.5, J = 2.2, 
77 = 1, at an accelerator age of 10 3 (filled circles), 10 4 years (open squares), and 10 s years (stars). Error 
bars are set to 10% of the count number, to mimic the expected error on the effective area for the array 
used in the counts determination. 

to continuous accelerators will present steeper profiles due to the freshly accelerated 
particles in the center of the source. Faster diffusion also leads to a larger extension. 
The profiles in Fig. 2 are normalized to their respective maximum counts, with flux 
decreasing with age in the impulsive case and opposite behavior in the continuous case. 

Acknowledgements We acknowledge discussions with Ana Y. Rodriguez-Marrero and many col- 
leagues of the CTA collaboration, especially Sabrina Casanova, Elsa de Cea del Pozo, Daniela Hadasch, 
Matthieu Renaud, Jim Hinton, Konrad Bernloehr, and Abelardo Moralejo. We acknowledge support 
from the Ministry of Science and the Generalitat de Catalunya, through the grants AYA2009-07391 and 
SGR2009-81 1, as well as by ASPERA-EU through grant EUI-2009-04072. 

REFERENCES 

1 . Albert J. et al. 2007 ApJ, 664L, 87A 

2. Aharonian F. et al. 20008 A&A, 481, 401 

3. Ackermann, M. et al. 2012 ApJ, 755, 22A 

4. Montmerle, T. 1979 ApJ, 231, 95M 

5. Aharonian & Atoyan 1996 ApJ, 309, 917 

6. Gabici et al. 2007 AsSS, 309, 365 

7. Rodriguez Marrero A. Y. et al., 2008, ApJ, 689, 213 

8 . Pedaletti et al. , 20 1 2, intended for A&A 

9. Ginzburg V. L., Syrovatskii S. I., 1964, The Origin of Cosmic Rays. Pergamon Press, London 

10. Berezinsky, V. S. et al.: Astrophysics of Cosmic Rays, North Holland (1990) 

11. Tarn, P. H. T., Wagner, S. J., Tibolla, O. & Chaves, R. C. G. 2010 A&A, 518A, 8T 

12. Actis, M. et al. (CTA Consortium) 201 1 ExA, 32, 193A 

13. Atwood, W. B., et al. 2009, ApJ, 697, 1071 

14. Goodman J. A., for the HAWC Collaboration 2010 ASPC, 426, 19G 

15. Bernlohr, K. et al. (CTA Consortium) 2012 Astroparticle Physics CTA Special Issue, in press 



